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BURNING OF SOLID ALUMINUM AND MAGNESIUM SPHERES IN 
HIGH-TEMPERATURE AND HIGH-VaOCITY GASES 
by Robert D. Ingebo 
Lewis Research Center 

SUMMARY 

Burning of solid aluminum and magnesium spheres in high-temperature and high- 
velocity gases, produced by the combustion of hydrogen and oxygen, was investigated 
experimentally. Metal spheres were 0.63 centimeter in diameter, and the flow environ- 
ment was varied to produce a Reynolds number of the sphere that varied from approxi- 
mately 6000 to 31 000. Burning or erosion started at the front surface of the metal 
sphere. The entire burning process was divided into two phases; the initial phase, dur- 
ing which the sphere was heated to the melting temperature of the metal with no erosion, 
was called a heating period. The heating period was followed by erosion with combustion 
stabilized in the wake of the sphere. 

The experimental heating and erosion times were correlated with Nusselt numbers 
by an equation derived from a model using a heat balance at the surface of the solid. The 
heat balance equation was based on the Reynolds and Nusselt numbers of the initial sphere 
because it was observed that the erosion rate was constant as the mass of the metal or 
effective diameter decreased. 


INTRODUCTION 

Rapid and efficient burning of high-energy fuels such as aluminum and magnesium is 
difficult to achieve. This is primarily due to combustion- inhibiting processes inherent 
in metal-fuel combustion. However, metal fuels still appear promising as a source of 
power. They have been used primarily as additives in solid-fueled rocket engines to 
improve the specific impulse and stability. 



Considerable effort has been made to determine the combustion mechanism for 
metals in quiescent oxidizing gases (refs. 1 to 3). Also, data have been obtained on the 
combustion of slurry fuels in jet engines (ref. 4). However, very little data is available 
on metal burni'r^ rates under actual combustor conditions and the influence of metals on 
combustion stability. The development of metal-fueled combustors requires new tech- 
niques of injecting, rapidly igniting, and burning metal particles to efficiently control the 
burning process. 

In this investigation of metal fuels, the main purpose was to determine mass and 
temperature histories of relatively large (0.6 3- cm diam) aluminum and magnesium 
spheres burning under simulated rocket combustor conditions. Such data are needed to 
derive a combustion model that could describe the actual burning process for metal fuels 
in a rocket engine and possibly explain the Influence of metal additives in controlling com- 
bustion instability. 


SYMBOLS 

Cp specific heat, g-cal/(g)(°C); J/(g)(^C) 

Dj. nozzle throat diameter, cm 

d sphere diameter, cm 

dm/dt erosion rate, g/sec 

latent heat of fusion, g-cal/g; J/g 

h heat-transfer coefficient, g-cal/(sec)(cm^)(°C); J/(sec)(cm^)(°C) 

k thermal conductivity, g-cal/(sec)(cm^)(°C/cm); J/(sec)(cm^)(°C/cm) 

m sphere weight, g 

Nu Nusselt number, hd/k 

o 

o/f mass mixture ratio of oxidant to fuel 

P pressure, psia; N/m abs 

Pr Prandtl number , M gCp ^ g/kg 

q heat-transfer rate, g-cal/sec; J/sec 

Re Reynolds number based on sphere diameter, dVgPg/jUg 

r sphere radius, cm 

T temperature , °C 

t time, sec 
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V 


velocity, cm/sec 
o!j absorptivity 

QI 2 emissivity 

2 

[X fluid viscosity, P; (N)(sec)/m 

Q 

p fluid density, g/ cm 

a Stefan- Boltzmann constant 

Subscripts: 

avg average for sphere lifetime 

C sphere center 

c convection 

e erosion 

fs front surface 

g gas 

i initial 

m metal 

mp melting point 

r radiation 

s surface 

t total 

Superscripts: 

— average 

' approximate 

APPARATUS AND PROCEDURE 

A rocket combustor, shown in figure 1, was used to determine burning histories of 
solid 0.63-centimeter-diameter aluminum and magnesium spheres inserted into high- 
velocity and high-temperature gases. The gas stream was produced by burning gaseous 
hydrogen and liquid oxygen injected 46 centimeters upstream of the spheres at mass mix- 
ture ratios of oxidant to fuel o/f of 1.7, 5.3, 8, 10, and 12. Three different exhaust 
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Figure 1. - Diagram of rocket combustor and auxiiiary equipment. 


nozzles were used having throat diameters of 2.46, 2.01, and 1.49 centimeters, which 
gave a combustion gas-stream velocity range of 6600 to 25 400 centimeters per second. 
These operating conditions resulted in a Reynolds number for the metal sphere ranging 
between 6000 and 31 000. At a given o/f condition, the chamber pressure was approxi- 
mately the same (within ±5 percent) for the three exhaust nozzle test conditions. Chamber 

5 5 2 

pressures varied between 13.6 and 18.4 atmospheres (1. 37x10 and 1. 86x10 N/m ) with 
different mixture ratios. 
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Spheres of aluminum (type 6061, 98 percent pure) and magnesium (type AZ31B, 

97 percent pure) were machined from 1/4-inch (0.635-cm) rods. A hole was drilled and 
tapped to the center of the sphere for mounting on a 1/16-inch (0. 159-cm) threaded 
stainless -steel rod. The spheres were mounted on the ceramic- coated tip of the probe, 
as shown in figure 1, and rapidly inserted into the center of the test section after equili- 
brium conditions had been obtained in the combustor. The tests were terminated by 
rapidly withdrawing the sphere back into the sheltered region cooled by nitrogen gas flow. 
The time measurement for each test began when the sphere passed out through the opening 
in the combustor wall and ended when it passed back through the same opening. The time 
required for the actuator to move the sphere 5. 1 centimeters was 2 milliseconds for either 
insertion or withdrawal. 

An iron-constantan thermocouple (0. 0127-cm-diam wire) was mounted in the center 
of the sphere to determine its temperature history at that point (r = 0). In obtaining mass 
histories, the quantity of metal eroded from the sphere was determined by weighing the 
sphere before and after each burning test. Attempts to determine the weight loss of a 
sphere by photographing it as it burned inside of the combustor proved unsuccessful since 
erosion was very nonuniform over its surface. 


RESULTS AND DISCUSSION 

A typical set of data showing mass and temperature histories of an aluminum sphere 
is plotted in figure 2(a). Actually, seven spheres weighing 0.350±0.0 l 3 gram were tested 
under identical conditions to obtain the complete history of a single sphere. One of the 
spheres was instrumented with a thermocouple junction installed in the center of the 
sphere to obtain the temperature history T^ as a function of time t. In the case of 
aluminum (fig. 2(a)), the temperature at the center rose from the initial temperature of 
20° C to the melting point (660° C) in approximately 250 milliseconds and then remained 
constant. In the case of magnesium (fig. 2(b)), the temperature at the center of the 
sphere rose to a constant value considerably below the melting point of magnesium in ap- 
proximately 150 milliseconds. 

Figure 2(a) also shows the mass history of an aluminum sphere. At approximately 
170 milliseconds after the sphere was inserted into the gas stream, it was observed, both 
photographically during the run and visually by inspecting the withdrawn sphere, that 
erosion was just beginning to occur at the front surface. Erosion at this surface is be- 
lieved to be due to the fact that the heat-transfer coefficient was considerably greater at 
the front surface than at the rear surface of the sphere. This is shown by the plots of 
calculated values of the front surface temperature Tfs against time t and the average 
surface temperature Tg against time t in figure 2(a). Thus, in 170 milliseconds, the 
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Data 
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n Sphere weight 

O Temperatu re at sphere center 

Front surface temperature 1 Calculated 

— ^ — Average surface temperature > from eq. (1) 

Temperature at sphere center J 

Experimental 

Solid symbols denote photographic data 



Time, t, msec 
(b) Magnesium. 

Figure 2. ~ Weight and temperature histories of 0.63-centimeter-diameter aluminum and 
magnesium spheres. Mass mixture ratio of oxidant to fuel, 1.7; nozzle throat diameter, 
2.01 centimeters. 
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front surface temperature reached the melting point (660^ C), whereas the calculated 
average surface temperature was 30° below the melting point. 

Spheres withdrawn at 205, 255, and 305 milliseconds are shown in figure 2(a). These 
photographs indicate that erosion was considerably greater at the front than at the rear of 
the sphere because of the higher heat-transfer coefficient at the leading surface. As a 
result, the metal spheres did not remain spherical during the erosion period. The mass 
history of the sphere showed that shortly after erosion began (called ignition) its weight 
decreased linearly with time. 


Heating of Metal Spheres 

Data on the sphere heating rate and the time required for it to reach the melting 
point at the front surface were analyzed in the following manner. 

Heating rates. - The heating rate was calculated by writing a heat balance at the sur- 
face of the metal sphere as 


2, dT 

Q = Q = jrr k — 
^m ^g,c m 


r=r. 


( 1 ) 


where q is the heating rate of the sphere and q^ is the rate of heat transfer by con- 
vection from the gas stream. Heat transferred by radiation was negligible (approx 
2 percent of q^ ^). As discussed in appendix A, the change in sphere temperature with 
time was obtained from equation (1) by simultaneously solving equations (Al) and (A3) at 
the desired radii (r = 0 and r = r ). 

Values of the heat-transfer coefficient were assumed to calculate several tempera- 
ture histories of the sphere. The best results shown in figure 2 were obtained with a 
heat-transfer coefficient h of 0.202 gram-calorie per second per square centimeter per 
°C (0. 847 J/(sec)(cm^)(°C)). Both figures show that agreement was good between calcu- 
lated and experimental values of the temperature at the center of the sphere T^. Also 
shown in figure 2 are values of T_, the average surface temperature of the sphere. 

Values of the average heat-transfer coefficient h are given in table I for a range of 
o/f conditions and corresponding values of the average Nusselt number, the Prandtl num- 
ber, and the Reynolds number. The aluminum and magnesium spheres were tested at 
five different o/f conditions. Combustion gas properties used in these calculations are 
discussed in appendix B. Gas velocities used for these tests are given in table H for a 
nozzle throat diameter of 2.01 centimeter. 
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TABLE I. - CALCULATED VALUES FOR HEAT-TRANSFER RATES 


llllllllll 


[Nozzle throat diameter, 2.01 cm.] 


Mass 
mixture 
ratio of 
oxidant 
to fuel, 
o/f 

Temper- 
ature of 
gas, 

^g’ 

Average heat-transfer coefficient, h 

Average 

Nusselt 

number, 

Nu 

Prandtl 
number , 
Pr 

Reynolds 

number. 

Re 

g- cal/(se c) ( cm C) 

J/(sec){cm2)(®C) 

1,7 

1545 

0.2020 

0.848 

132. 5 

0.526 

21 000 

5.3 

2977 

.1295 

.543 

93.0 

.595 

12 130 

8 

3198 

.0967 

.406 

94.0 

.758 

10 850 

10 

3144 

.0904 

.379 

98.5 

,744 

10 400 

12 

3054 

.0812 

.341 

94.5 

.734 

10 380 


TABLE II. - EXPERIMENTAL RESULTS 


Nozzle 

Mass 

Velocity 

Reynolds 

Average 


Magnesium 



Aluminum 


throat 

diameter, 

Df 

cm 

mixture 
ratio of 
oxidant 
to fuel, 
o/f 

of gas, 

y 

cm/sec 

number, 

Re 

Nusselt 

number, 

Nu 

Heating 

time, 

‘s> 

msec 

Average 

surface 

temper- 

ature, 

Ts> 

°C 

Average erosion 
rate, 

(dm/dt)^^g, 

g/sec 

Heating 

time, 

‘s- 

msec 

Average 

surface 

temper- 

ature, 

Average erosion 
rate, 

(dm/dt)^^g, 

g/sec 







Experi- 

mental 

Calcu- 

lated 


^S’ 

°C 

Experi- 

mental 

Calcu- 

lated 









(eq. (10)) 




(eq. (10)) 

1,49 

1.7 

9 270 

11 500 

90,5 

175 

555 

1.48 

1.42 


— 


— 


5.3 

8 150 

6 650 

67.7 

115 

540 

2. 50 

2.64 



— 

— 


8 

7 530 

5 950 

67.3 

140 

524 

1.84 

1.78 

--- 

___ 

— 

— 


10 

6 830 

5 700 

65. 5 

160 

497 

1.54 

1.44 



— 

— 


12 

6 580 

5 680 

65.3 

170 

479 

1.32 

1.26 

— 

-- 

— 

— 

2.01 

1.7 

16 900 

21 000 

132.5 

110 

525 

1.84 

2.02 

170 

630 

2.07 

2.29 


5.3 

14 850 

12 130 

93.0 

80 

526 

2.80 

3.05 

120 

620 

3.27 

3.63 


8 

13 730 

10 850 

94.0 

90 

493 

2.15 

2.34 

150 

610 

2.89 

2.91 


10 

12 450 

10 400 

98.5 

95 

476 

2.00 

2.06 

160 

603 

2.71 

2.61 


12 

12 000 

10 380 

98.3 

100 

450 

1.77 

1.78 

175 

595 

2.43 

2.31 

2.46 

1.7 

25 400 

31 600 

169. 5 

90 

520 

2.23 

2.56 

135 

625 

3. 12 

2.94 


5.3 

22 300 

18 200 

123.5 

60 

497 

3.54 

3.83 

95 

614 

4.38 

4. 85 


8 

20 600 

16 300 

125.0 

70 

481 

3.08 

3.03 

115 

600 

3.76 

3.88 


10 

18 700 

15 600 

118.5 

75 

447 

2.38 

2.65 

125 

592 

3.36 

3. 15 


12 

18 000 

15 570 

118.5 

80 

431 

2.30 

2.08 

135 

585 

3.07 

2.70 


8 




mil II I 


I 


The heat-transfer coefficients were correlated with the dimensionless groups as fol- 
lows: McAdams (ref. 5) relates the heat-transfer coefficient for solid spheres at large 
Reynolds numbers with the Nusselt number and Reynolds number as 


— = Nu = 0.37 


( 2 ) 


However, the theoretical minimum value for the Nusselt number is 2 (when Re =0). 

0 3 

Heat-transfer relations have also generally shown that Nu Pr ' . Therefore, the pre- 
ceding expression was modified to 


Nu = 2 + CPr'^-^Re®'® (3) 

where C is a proportionality constant. To evaluate C, data in table I were plotted as 
shown in figure 3, and the following expression was obtained for solid spheres: 

Nu = 2 + 0.4 Pr^-^Re®-® (4) 

Properties used in equations (2) to (4) are evaluated at the mean temperature of the com- 
bustion gases and surface and for static pressures. 


Mass mixture 
ratio of oxidant 



Figure 3. ’ Correlation of Nusselt number with Prandtl 
and Reynolds numbers. 
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Aluminum spheres 

1 


Mass mixture Reynolds 
ratio of oxida nt n u mbe r, 
to fuel, Re 
o/f 

O 1.7 31600 

□ 5.3 18 200 

O 8 16 300 

> 10 15 600 

F 12 15 570 


Magnesium spheres 


\ ^ Calculated total time 

n Experimental 

\\ ^ 'A V Calculated from eqs; (8) and Ul) 




Solid symbols denote photographic data 




40 80 120 160 200 240 280 

(a) Nozzle throat diameter, 2.46 centimeters. 


Aluminum spheres 


° 21000 

\\ K\\\ \\ □ 5.3 12130 

\v \\^ V O 8 10 850 

lA \ \ [> 10 10 400 

V \\ \ 'a ^ Calculated total time 

\\ \\ \ Experimental 

Magnesium spheres U U ^ \\ Calculated from eqs. (8) and (11) 

\ Solid symbols denote photographic data 


Mass mixture Reynolds 
ratio of oxidant number, 
to fuel, Re 
o/f 












100 140 180 220 260 300 "" 340 380 

Time, t, msec 

(b) Nozzle throat diameter, 2.01 centimeters. 

Figure 4.^- Erosion histories. 


Ml i i mi i Minniii || IHIIH || 


Heating times. - From high-speed motion pictures of the spheres, it was observed 
that melting first occurred at the front surface. Thus, the heating time t^ was deter- 
mined photographically from the time of insertion of the sphere into the gas stream to the 
time that melting was first observed at the front surface. Experimental values of the 

heating time t are recorded in table n and plotted as the solid symbols in figure 4 at the 
s 

breakpoint in the mass history of the sphere. 



Mass mixture 

Reynolds 


ratio of 

number, 


oxidant to fuel, 
o/f 

Re 

o 

1.7 

11500 

□ 

5.3 

6 650 

O 

8 

5 950 

D> 

10 

5 700 

£7 

12 

5 680 

O 

Calculated total time 



Figure 4. - Concluded. 


In order to calculate the heating time t from the heating rate (eq. (l) or (Al) and 
(A3)) , it is necessary to know the average surface temperature Tg at the time t^ (when 
Tfs = Tjj^p). The average surface temperature when melting first occurred varied with 
flow conditions. Therefore, it was assumed that the nondimens ional temperature ratio 
T /T _ was a function of the average Nusselt number. Also, T was considerably lower 

S III p o 
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for magnesium than it was for aluminum, at the same value of Nu. Consequently, it was 
also assumed that T„/T was a function of the ratio of gas to liquid-metal density 

o IXip 

Pg/Pm" The fallowing expression was therefore assumed for calculating Tg at the time 
^ g m 


^ = 1 + amH ( 5 ) 

Ts \Pm) 

This is an empirical correlation that might not be valid for other metals. 

The exponent m was evaluated by plotting (T^ /T„) - 1 against Nu (as shown in 

O _ -1 n 

fig. 5), and m was found to be 0.7. Similarly, a plot of (Tj^^p/Tg )Nu against Pg/Pm 


6r" 





Mass mixture 
ratio of oxidant 
to fuel, 
o/f 

O 1.7 

□ 5.3 



1 I I . J . 1 I- 1 I 

40 60 80 100 200 

Average Nusselt number, Nu 

Figure 5. - Determination of exponent m for average Nusselt 
number. 


(shown in fig. 6), gave an n of 3. Finally, T^p/Tg was plotted against^ Nu^* "^(Pg/p^^) 3 
(fig. 7 ) , and the proportionality constant C was determined to be 1. 1x10 . Thus, equa- 
tion (5) was rewritten as 


12 





’^ = 1 ^1.1X10^ 


which gives the average surface temperature Tg when the leading edge reaches the melt- 
ing temperature, defined as the end of the heating time t . This is the time at which 

s 

heating ends and erosion begins . 

Although machine calculations using equations (1) and (6) were used to determine t 
— s 

at Tg, a simpler method was derived from the following heat balance at the surface of the 
sphere; 


/T^\ dT 

»dik^(Tg-T3)NUj=mCp_ 


which may be integrated to give 


,3 Tg - 

irdjkgNUjT^ T - T3 


TABLE III. - COMPARISON BETWEEN CALCULATED AND EXPERIMENTAL 
RESULTS FOR MAGNESIUM AND ALUMINUM SPHERES 
[Nozzle throat diameter, 2.01 cm.] 


Metal Mass 
mixture 


Magnesium 1.7 
5.3 


Aluminum 



ime, msec 

Erosion time, t^, msec 

Total tim 

Experimental, 

exp 

Equation (11) 

Experimental 

Calculated 

110 

114 

125 

225 

80 

76 

83 

156 

90 

99 

108 

191 

95 

112 

115 

207 

100 

130 

130 

228 

180 

153 

174 

327 

120 

97 

120 

222 

150 

120 

138 

270 

160 

134 

150 

295 

175 

151 

169 

322 
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where t' is the approximate heating time and T. is the initial sphere temperature. 

o 1 

The agreement between values of t’ , calculated from equation (8), and experimental 

s 

values of t was good, as shown in table III. Thus, t’ is the approximate heating time 
s s 

required for the average surface temperature to reach T and erosion to begin. 

s 

Erosion of Metal Spheres 

Motion pictures and eroded spheres indicated that molten metal was stripped from the 
surface of the sphere by the high-velocity gas stream and that the atomized liquid was va- 
porized and ignited, and combustion was stabilized in the wake of the sphere. A model 
for this erosion process is shown in figure 8. 




Figure 8. - Combustion model for eroding metal spheres. 


In this portion of the study, the main objective was to determine erosion (mass re- 
moval) rates at the surface of the sphere and to compare them with calculated values. 
Calculations were made based on a heat balance and mass balance at the surface of the 
sphere. The use of a mass balance proved unsuccessful since the mass removal or ero- 
sion rate was not limited by the diffusion rate of oxidant gases to the surface of the sphere. 
Thus, a heat balance at the surface of the sphere was used to analyze the erosion data. 

Erosion rates. - Applying a heat balance at the surface to the erosion of a metal 
sphere in high-temperature and high-velocity gas streams gives 

qL = "ig 

where qj^ is the rate at which heat is gained by the eroding liquid and q^ is the rate of 
heat transfer from the gas stream. In deriving this expression, the heating rate of the 
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metal was assumed to be negligible since it decreased rapidly when erosion began, and 
heat transfer from the flame was assumed negligible. Figure 4 shows that the average 
erosion rate remained approximately constant with time. This result was 

attributed to the fact that as the metal eroded it did not remain spherical in shape. Thus, 
Nu could not be accurately determined during erosion. In fact, the curvature of the front 
surface became less convex, which tended to increase the front surface Nusselt number 
although the sphere diameter decreased with time. Therefore, the initial Nusselt number 
before erosion started was used to define the heat transfer to the surface during erosion. 
Because of the nonuniformity of the surface temperature of the sphere, the average value 

of the surface temperature T„ was used to evaluate the sensible heat gained by the metal 

s 

as it was melted and eroded by the gas stream. Thus, T was also used to evaluate the 

s 

temperature drop across the gas film surrounding the sphere. A heat balance at the sur- 
face of the sphere was written as 




+ c (T 
p,m' mp 



7^d,k„(T„ - T„)Nu. 
i g'‘ g s^ 1 


( 10 ) 


where (dm/dt) is the average erosion rate, is the heat of fusion of the metal, and 
d^ and Nuj^ are the diameter and average Nusselt number, respectively, based on the 
initial sphere diameter. Values of calculated from equation (10) and 

plotted as shown in figure 4. 

Erosion times. - The erosion time t^ may be obtained by integrating equation (10): 


t 


e 


pi<[- 




6k (T - T, 
g' s s 


mp "^s^] 

)W.“ " 


( 11 ) 


Values of t^ calculated from this equation agree reasonably well with experimental 
values of t^, as shown in table m. 

Total heating and erosion time. - The total time required for heating and erosion 
of metal spheres can be calculated from 




( 12 ) 


By substituting equations (8) and (11) into equation (12), the following expression was 
obtained: 
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(13) 


t- 

i: — 

t XT., 


6k Nu. 
g 1 


T /T 

mp In I g 


kT - T 
^ g .sy 


^ '‘' ^p, m^^ mp ~ 


^p,m^^g ■ 


A comparison of values of calculated from equation (13) and experimental values is 
shown in figure 4 and table IIL Although no direct comparison could be made with the 
data given by Macek in reference 6, it could be assumed that Nu — 2 for the 0.0032 centi- 
meter diameter aluminum sphere. This gives t^ = 1.26 milliseconds, which agrees fairly 
well with his data. 


Burning Process with Eroding Metal Spheres 

High-speed motion pictures of the spheres showed that soon after erosion was initi- 
ated the molten metal ignited and combustion was stabilized in the downstream wake of the 
sphere. A combustion model illustrating this burning with solid spheres is shown in fig- 
ure 8. As this burning model indicates, there was no evidence of solid-phase combustion 
at the surface of the sphere, and no oxide shell was formed that would inhibit erosion or 
burning. The eroded molten metal appeared to be vaporizing and combustion was stabil- 
ized in the wake of the sphere with no appreciable amount of heat being transferred from 
the flame to the surface of the sphere. The luminosity of the flame was greatly increased 
when the erosion rate was increased and also when the o/f was increased. 


SUMMARY OF RESULTS 

The burning of solid aluminum and magnesium spheres (0.63-cm diam) in high- 
temperature and high-velocity gases, produced by the combustion of hydrogen and oxygen, 
has shown the following; 

1 . The front surface of a metal sphere was heated to its melting point before it began 
to erode and burn. There was no evidence of solid-phase burning at the surface of the 
sphere, and therefore no oxide shell was formed to inhibit the erosion process. Metal 
combustion was stabilized in the wake of the sphere and did not appreciably affect the 
erosion rate by transferring heat to the surface of the sphere. 

2. Experimental heating and erosion times agreed well with values calculated from a 
heat balance at the surface of the sphere. The heating time for a metal sphere was 
approximately one-half as long as the total time required for complete erosion of the 
sphere. 

3. The total burning time t^. (heating time plus erosion time) was correlated with the 
average Nusselt number Nu: 
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2 r 

w, /T„ - TA 

6kgNu [t^ \Ts-V 

where is the fluid density of the metal, d, is the initial sphere diameter, c „ is 

^ ^ PjHi 

the specific heat of the metal, is the thermal conductivity of the gas, is the tem- 
perature of the metal, Tg is the average surface temperature of the sphere, Tg is the 
temperature of the gas, is the initial temperature of the metal, is the latent heat 
of fusion, and is the melting point of the metal. 

4. The average surface temperature of the sphere Tg was correlated with the aver- 
age Nusselt number Nu and the ratio of gas to liquid-metal density p^/p^'- 

3 

= 1.0 + 1. 10X10^ 

^s \Pm / 

5. The following correlation of the average Nusselt number with the Prandtl and 
Reynolds numbers was obtained: 

Nu = 2 + 0.4Pr°-^Re°-® 

This agreed well with data in the literature for heat transfer from air streams to solid 
spheres. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 11, 1970, 

128-31. 


^f ^p,m^^mp ~ ^s^ 
Vm<Tg-T,) , 
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APPENDIX A 


HEATING OF SOLID METAL SPHERES IN GAS STREAMS 


The conduction of heat in a solid metal sphere may be expressed in spherical coor- 
dinates as 


J.i./,2 3T\ jT ,Ai 

where k^/p^^Cp is the thermal diffusivity or ratio of the thermal conductivity to the 
product of the density and specific heat of the sphere and r is the radial distance of a 
point from the center of the sphere at temperature T and time t. 

The heat balance at the surface of a metal sphere suspended in a high-temperature 
and high-velocity gas stream may be written as 


V ^g,c ^g,r 


(A2) 


where q is the heating rate of the sphere and q and q are rates of heat trans- 
fer from the gas stream by convection and radiation, respectively. Since 


q 


m 


= 77r2k^ 


dT 

dr 




and 


= irr 




T„ - Q!oT 


2"S; 


equation (A2) becomes 


dT 

dr 


= (T„ - T ) 

k g s' 
r=rg m 




(A3) 
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where h is the heat-transfer coefficient at the surface of the sphere, a is the Stefan- 
Boltzmann constant, and and oig are the absorptivity and emissivity of the' polished 
metal, respectively. 
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APPENDIX B 

COMBUSTION GAS PROPERTIES 

Combustion gas properties were calculated for equilibrium conditions in the combus- 
tor by the method given in reference 6. The results are given in table IV for the three 
exhaust nozzles used in this study. 

The viscosity and thermal conductivity of the combustion gas stream were also deter- 
mined by this method for a range of gas temperatures. The results of these calculations 
are plotted in figure 9. 


TABLE IV. - COMBUSTION GAS STREAM PROPERTIES 


Mass mixture 
ratio, 
o/f 

Gas pressure, 

Gas temperature, 
^g’ 

Gas density, 
^g’ 

psia 

N/m^ abs 

o 

o 

g/cm^ 

1.7 

270 

18,6X10^ 

1818 

6.70 

5.3 

255 

17.5 

3250 

7,92 

8 

230 

15.8 

3471 

8.58 

10 

215 

14.8 

3417 

9.18 

12 

200 

13.7 

3327 

9.57 


4r— 


e 


2h- 



10 


-5 



Figure 9. - Combustion gas properties calculated by method given in reference 7. 
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